Light output response of the LVD liquid scintillator to neutron-induced nuclear recoils 
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The organic liquid scintillator used in the LVD experiment (INFN Gran Sasso National Labo- 
ratory) has been exposed to an Am-Be neutron source to measure the light response function for 
neutron energies in the region from about 4 to 11 MeV. 

A full Monte Carlo simulation, incorporating the detector response, is used to generate neutron 
scattering spectra which are matched to the observed ones to determine the quenching factors. 
The obtained light output response is well described by the semi-empirical Birks model. 
The results, consistent with those obtained by other authors using similar hydrocarbonate scintilla- 
tors, can be of interest for fast and high-energy neutron spectroscopy that could be performed with 
this detector. 
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INTRODUCTION 



The Large Volume Detector (LVD), in the INFN Gran 
Sasso National Laboratory (Italy) , below an average 3600 
m water equivalent rock overburden, is a 1 kt liquid scin- 
tillator detector whose major purpose is monitoring the 
Galaxy to study neutrino bursts from gravitational stellar 
collapses [1]. LVD consists of an array of 840 scintillator 
counters, 1.5 m 3 each, divided in three, identical and in- 
dependent towers. Each counter is viewed from the top 
by three 15 cm photomultiplicr tubes (PMTs) [2]. 
The LVD liquid scintillator is a mixture of aliphatic 
and aromatic hydrocarbons (Cnld^n with n = 9.6) also 
known as White Spirit. For historical reasons 73% of the 
LVD scintillator was produced in Russia and has an aro- 
matic hydrocarbons content of ~16%, while the remain- 
ing amount has been produced in Italy and contains ^8% 
of aromatics. The scintillator with the higher percentage 
of aromatics presents a light yield of about 20% higher. 
We do not expect this slight difference on the scintillator 
composition to have an impact on the quenching factors 
therefore we decided to measure the proton light output 
response in the liquid scintillator produced in Russia that 
is present in the majority of counters (LVD-scint). 
In addition, since in 2005 one of the LVD counters (with 
aromatic hydrocarbons content of ~8%) was doped with 
gadolinium [3] to improve the signal to noise ratio in neu- 
tron detection, we decided to repeat the measurement 
also on this particular scintillator (LVD-Gd). 
In organic scintillators the light output response, the en- 
ergy emitted as fluorescence, is non-linear with respect to 
the energy deposited by ionizing particles (this is partic- 
ularly true for highly ionizing particles). This behavior 
is attributed to quenching of the primary excitation by 
the high ionization density along the particle track which 
leads to a decrease of the scintillation efficiency. In the 
absence of quenching the light yield is proportional to 
the energy released, which may be written in differential 



form as: 



dL _ dE 
dx dx 



(1) 



where ^ is the fluorescent energy emitted per unit path 
length, S is the scintillation efficiency and 4^ is the spe- 
cific stopping power. To describe the quenching effect, 
Birks [4] proposed a semi-empirical relation in which the 
differential light output is reduced for high energy loss: 
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where the product kB is usually treated as a single pa- 
rameter and it is known as Birks factor. 
The quenching factor for ions is defined as the ratio of 
light yield of ions to that of electrons of the same energy. 
Integrating eq. 2 it can be written as: 
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EXPERIMENTAL SETUP AND PROCEDURE 

Two scintillator detectors of the same dimensions of 
the counters of the LVD experiment, have been used to 
perform the measurements. Each detector consists of a 
1 x 1 x 1.5 m 3 stainless steel tank filled with 1.2 ton 
of liquid scintillator and viewed from the top by three 
PMTs PHOTONIS XP3550 of 5" diameter. The sig- 
nals of each PMT are recorded by a CAEN VI 724 wave- 
form digitizer, with 100 MHz of sampling frequency and 
137 /j,V of resolution per bit. The board hosts a Field 
Programmable Gate Array (FPGA) for each input chan- 
nel. The data stream is continuously written in a circu- 
lar memory buffer, which is frozen by the FPGA when 
the trigger occurs. The board implement also the Zero 
Length Encoding (ZLE) algorithm for data reduction. 



This algorithm allows to set a threshold, which could be 
different from the trigger threshold, and discard the data 
under threshold saving time in both data transfer and 
processing. The acquisition window length, the num- 
ber of samples pre- and post-trigger, the threshold and 
the level of coincidence among different channels are pro- 
grammable features. 

In the center of the detector we have inserted an 
241 Am 9 Be source encapsulated in a waterproof housing 
and held in place by a stainless steel rod. 
The radionuclide 241 Am decays emitting an alpha parti- 
cle of 5.6 MeV. Under bombardment by alphas, 9 Be un- 
dergoes a nuclear reaction forming the excited compound 
nucleus 13 C*, which then decays mostly producing either 
12 C at ground state or at the first excited level [6]: 
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E n < 6.5 MeV 
E 7 =4.44 MeV (5) 



In the latter case the neutron is accompanied by a de- 
excitation gamma-ray of energy 4.44 MeV. The source 
used in these measurements emits about 10 neutrons s _1 . 
The emitted neutrons are thermalized by elastic scatter- 
ing on the hydrogen atoms of the organic scintillator. 
During the slowing down process, a fraction of the neu- 
tron kinetic energy is converted to photons which are 
detected by the PMTs giving the prompt signal. After 
being thermalized, neutrons are eventually captured by 
■"■H, na *Gd (when present) and with lower probability by 
12 C and 56 Fe nuclei producing the delayed signal. The 
mean neutron capture time depends on the target com- 
position (n-capture cross section) and on the detector ge- 
ometry because neutrons which migrate farther from the 
source can escape from the detector leading to a shorter 
mean capture time. 

ANALYSIS AND RESULTS 

In the present work we will refer to the light output 
in terms of the scintillator response to electrons, as this 
can be taken to be linear (within ~1% for electrons of 
energy E > 2 MeV). The measured signals will then be 
expressed in MeV electron recoils equivalent (MeV ee ). 
Assuming that our scintillator has a response function 
to nuclear recoils consistent with that of similar organic 
scintillators, we should expect signals up to ^2 McV ee 
for neutrons of 6.5 MeV which elastically scatter protons. 
Neutrons of lower energy are accompanied with a 4.44 
MeV gamma emission, according to equations 4 and 5. 
Neutron captures on H are followed by the emission 
of the 2.22 MeV 7 quantum from the deuterium de- 
excitation; while captures on Gd originate 7 cascades 
of total energy ^8 MeV. As we want to measure the 



prompt signals with a threshold as low as possible, we 
decided to trigger on the delayed pulse and then look 
backward in time to find the prompt one. The trigger 
condition to enable samples storage was chosen to be 
a 3-fold coincidence of the PMTs signal exceeding a 
threshold of 1.8 McV ee for the counter with LVD-scint 
and 3.2 MeV ee for the counter with LVD-Gd. The 
background rate at these thresholds, for a counter placed 
underground, is about 50% and 2% of the neutron rate 
from the source. 

Neutron events are selected requiring that the delayed 
pulse must be accompanied, in a time window of 320 
/l*s preceding the trigger, by a second signal (prompt) 
detected with a ZLE threshold of 0.8 MeV ee for both 
the detectors. This additional condition helps to get rid 
of the background enanching the signal to noise ratio of 
about a factor ten. 

The highest-energy prompt signals are large enough to 
trigger by themselves, in this case the delayed pulse will 
be left out of the acquisition window. To avoid this 
situation the acquisition window has been increased to 
640 /us moving the trigger in the middle of it. In this 
way we can consider the trigger pulse as due to nuclear 
recoil or nuclear capture, depending on the presence 
of a pulse in the following or in the preceding 320 /is 
respectively. In turn the pulse in the pre-trigger region 
can be regarded as a recoil signal and used to extend 
the recoil energy spectrum down to 0.8 MeV ee . 
In fig. 1 we show the distribution of the absolute value 
of the difference between the trigger pulse and the 
pre-trigger or post-trigger pulse if present in the same 
event, for both the scintillators measured. Fitting the 
distribution with an exponential function we found, for 
LVD-scint, a mean neutron capture time of 210 (is, well 
in agreement with previous measurements taking into 
account that the neutron source was placed in the center 
of the detector. Concerning the counter filled up with 
LVD-Gd, doped at a Gd concentration of 0.93 g/1, we 
found a mean neutron capture time of 27/is. Also in this 
case the measured value is in agreement with previous 
measurements and simulations [3]. In the same figure, 
gray points represent the background level, obtained 
applying the same analysis to a dataset collected after 
removing the Am-Be source. 

We developed a Monte Carlo simulation based on the 
GEANT4 toolkit [5] implementing the detector descrip- 
tion and the Am-Be source energy spectrum taken from 
[12]. 

GEANT4 has high precision transport models 
(G4NeutronHPElastic and G4NeutronHPInelastic) 
available for neutrons of energies up to 20 MeV and for 
all materials. The high precision models are based on the 
Evaluated Nuclear Data Files (ENDF/B-VI) [7], which 
contains tabulated cross-sections for elastic reactions and 
the different inelastic final states (e.g.: n7, np, nd, etc.). 
The high-quality ENDF/B-VI data generally allows for 
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FIG. 1: Distribution of the absolute value of the time dif- 
ference between the trigger and pre- or post-trigger in the 
same event (blu points) compared with the background level 
(in gray). Shapes are, as expected, flat for the background 
(fit function: y— kbck) and exponential for the data with the 
neutron source (fit function: y— Ae~ l ^ T + k^ck) ■ Top and bot- 
tom panel are relative to LVD-scint and LVD-Gd scintillator 
respectively. 



accurate simulations of these low-energy neutrons [8] . 
The outcome of the simulation is shown in fig. 2 where 
the energy spectra of the de-excitation 7-ray following 
neutron capture are compared with experimental ones. 
The 2.22 MeV gamma emitted after neutron capture 
on 1 H, clearly visible for both the detectors, has been 
used to determine the energy scale with an uncertainty 



FIG. 2: Energy spectrum of the neutron capture signals 
(background subtracted) (in black), compared with simula- 
tion (in red). Top and bottom panel are relative to LVD-scint 
and LVD-Gd scintillator respectively. 



of ±1.5%. 

The spectrum relative to the unloaded scintillator 
exhibits the peaks due to neutron captures on 12 C (4.95 
MeV), and 56 Fe (7.65 MeV), which are overwhelmed in 
the other detector by the presence of gadolinium. 
The light output produced by the nuclear recoils was 
evaluated from the eq. 3, taking into account the con- 
tribution of both the ions 1 H and 12 C, even if the latter 
was rather negligible when compared to the former. 
The stopping power dE/dx in liquid scintillator, from 
which the quenching factor depends, was taken from the 



ESTAR database [9] for electrons, whereas for protons 
and heavier ions it was calculated with the SRIM code 
[10]. 

For a specific material, the eq. 3 has only one free pa- 
rameter, kB, the so-called Birks factor, which has been 
determined by fitting, with the least-squares method, 
the simulated energy spectrum of nuclear recoils to 
the experimental one. In the fitting procedure, the 
comparison between experimental data and simulation 
was performed by re-weighting Monte Carlo events for 
different choice of the kB parameter. The parameter 
estimation and the goodness of fit are based on the x 2 
test, modified for comparing weighted and unweighted 
events [13]. 

The result of the fitting procedure is shown in fig. 3, 
where experimental and simulated date are represented 
with black and red line respectively. The two spectra 
diverge at E<0.8 MeV ee because of the detection thresh- 
old. The first peak at E~1.5 MeV ee is due to neutron 
elastic scattering, whereas the second peak is due to the 
4.44 MeV gamma from 12 C(n,n'7) in addition to elastic 
scattering of lower energy neutrons (E„ <6.5 MeV). 
The probability to get by chance a \ 2 higher than the 
observed is 34% for the Russian scintillator and 17% for 
the Italian one. 

The minimum \ 2 corresponds to a Birks param- 
eter value kB = 0.0125 g cm~ 2 MeV~ l and 
kB = 0.0140 g cwT^MeV' 1 for LVD-scint and 
LVD-Gd respectively. By letting \ 2 vary by one unit we 
found the standard error interval, which is ±1% in both 
cases. In addition to the statistical error we estimated a 
systematic contribution of ±5% coming mostly from the 
uncertainty in the energy calibration. 
According to the author of the work [11], the Birks 
factor is independent of the particle type and can then 
be used afterwords to calculate quenching factors for 
other particles and in other energy regions. 
In fig. 4 we show the resulting quenching curve for 
protons, alpha particles, and carbon ions. 



CONCLUSION 

Two different detectors of the LVD experiment have 
been exposed to an Am-Be (a,n) source, with the pur- 
pose of measuring the light output response to nuclear 
recoils. 

The delayed signal due to de-excitation 7-ray following 
neutron capture has been used to tag the events and ef- 
fectively reject the gamma background. 
The analysis is based on comparison of the experimental 
nuclear recoil spectra with Monte Carlo simulation. We 
found that the light output response is well described 
by the Birks model and a Birks factor kB = (0.0125 
± 0.0004) g cm^MeV' 1 and kB = (0.0140 ± 0.0007) 
g cm~ 2 MeV~ 1 was obtained from the two scintillators 
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FIG. 3: Energy spectrum of the prompt signals after back- 
ground subtraction (in black), compared with the simulation 
(in red). The detection threshold is at about 0.8 MeVee. 
The value of the Birks factor which minimize the x 2 i s 
kB = 0.0125 g cm' 2 MeV' 1 for LVD-scint (top panel) and 
kB = 0.0140 g cm^MeV' 1 for LVD-Gd (bottom panel). 



LVD-scint and LVD-Gd respectively. 
Quenching factors for ions in scintillator have then been 
calculated up to 100 MeV. 

This result can be of interest for fast and high-energy 
neutron spectroscopy performed with the LVD detector. 
It can be directly applied to 73% of the whole LVD ar- 
ray (Russian scintillator, i.e. LVD-scint) concerning the 
remaining amount of counters, filled up with Italian scin- 
tillator, we have not yet made a direct measurement of 
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FIG. 4: Quenching factors for protons, alpha particles and 
carbon ions obtained from the eq. 3 using the best fit param- 
eter kB = 0.0125 g cm^MeV' 1 for LVD-scint (solid lines) 
and kB = 0.0140 g cm^MeV' 1 for LVD-Gd (dashed lines). 



the quenching factors. Nevertheless, we could probably 
consider the results coming from the Italian Gd-doped 
scintillator as a lower bound for the undoped one. 

We are very grateful for helpful discussion with Dr. 
A.Ianni, Dr. M.Laubenstein and Dr. L.Pandola. 
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